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Abstract 



We investigate how Drell-Yan processes can be used to measure the magni- 



^ ■ tude of flavor symmetry violation in the proton sea. We examine the utility 



of the following beams: protons, charged pions, and charged kaons. In each 
case we present an approximate expression for the Drell-Yan asymmetry. Us- 
ing currently available parton distributions, we locate those kinematic regions 
which provide the greatest information on the quantity d^(x) — ^^{x). If suf- 
ficiently intense K^ beams were available, they could provide an efficient 
measurement of this quantity. Finally we present and discuss sets of sum 

rules for the Drell-Yan processes. 
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I. INTRODUCTION 

At the present time there is considerable interest in the flavor structure of the nucleon 
sea |jl| . This is largely a consequence of recent experimental results, such as the discovery 
by the New Muon Collaboration (NMC) of a violation of the Gottfried sum rule [Q. The 
most likely explanation of the NMC measurement was a significant flavor asymmetry in the 
proton sea, i.e. (i^(x) 7^ u^{x). Although there is no compelling reason why one should have 
(F{x) = vP{x)^ and in fact due to the Pauli principle one expects /q dx (F{x) — vP{x) > 0, 
phenomenological fits to parton distributions prior to the NMC measurement had typically 
set these sea distributions equal. 

Ellis and Stirling |^ pointed out that this asymmetry could be measured in proton- 
induced Drell-Yan (DY) processes on hydrogen and deuterium targets. They advocated 
measurements taken at center-of-mass rapidity y = 0, i.e. Xi = X2, where Xi {X2) is ap- 
proximately the momentum fraction carried by the projectile (target) quark or antiquark 
which annihilates in the DY process. These measurements were carried out by the NA51 
collaboration at CERN [ffl, which obtained the value 



u^{x) 



= 0.51±0M{stat)±0.05{syst). (1.1) 



dP{x 

The vP /dP distributions have recently been measured in the E866 experiment at FNAL 
0, through proton- induced DY processes on hydrogen and deuterium. The E866 group has 
measured uP jd^ over a range of X\ and X2 values. With the Main Injector upgrade scheduled 
for completion in early 1999 at FNAL, there is a revival of interest in the physics issues 
which could be addressed with the beams of kaons that could become available. We think 
that there is interesting physics which could be done in that program. 

In this paper we will review the prospects for measuring the flavor asymmetry in the 
nucleon sea through DY processes induced by protons or charged mesons [K^ and vr^). We 
will examine the sensitivity of various cross sections to the up/down antiquark ratio in the 
proton. We will also look at the sensitivity to the kinematic region for both projectile and 



target. We will discuss what measurements are needed if the DY process is to be used to 
pin down the largely unknown kaon structure functions. Finally we will present and discuss 
sets of sum rules for the DY processes. 

II. ANALYTIC APPROXIMATIONS FOR DRELL-YAN ASYMMETRY 

The general form of the DY cross section is: 

In Eq. ( p.l|) , a is the fine-structure constant, M^ = 5x1X2 is the dilepton mass squared, e^ 
is the charge of parton of flavor j, q and q are the quark and antiquark distributions, and 
the superscripts B and T refer to the beam and target, respectively. 

In order to probe differences between up and down antiquarks in the proton, one measures 
the ratio of DY cross sections on hydrogen and deuterium Qj. If one makes the impulse 
approximation, the DY cross section on deuterium {D) will just equal the sum of DY cross 
sections on the free proton {p) and neutron (n). For a beam particle B the asymmetry 
Ag^ = 2aBp/o-BD — 1 then becomes 

Aq (^Xi,X2j — jj£)Y( T"; tjDY ( V' l^"^J 

Because A^"^ is a ratio of cross sections, one expects some systematic errors to cancel. In 
subsequent sections, we will suppress for simplicity the dependence of A^^ on xi and X2. 

In comparing DY cross sections on protons and neutrons, we invoke charge symmetry 
for the nucleon parton distributions. This involves, for example, setting 

m"(x) = d^{x) and rf"(x) = mP(x), (2.3) 

and similarly for the antiquark distributions. We can therefore define all nucleon parton 
distributions in terms of those in the proton. Charge symmetry violating [CSV] effects have 



been estimated by Londergan, Thomas and collaborators PJ^ . For these processes the CSV 
effects should never be greater than 1% of the relevant amplitudes. 

Contributions to v4^^ which involve an up quark from the projectile will then be directly 
proportional to u'p{x) — dP^x). Since up quark terms will be four times as large as down 
quark terms, due to the square of the quark charge, such terms should make substantial 
contributions to the measured asymmetries. For this reason A^ should be an excellent 
place to measure up/down differences in the proton sea. This does not mean that all 
projectiles would be equally good probes of vF{x) — (F{x). In fact p and K^ turn out to be 
better than K^ , 7r+ and tt". Reasons for this are discussed in section IV. 

In subsequent sections, we make the following approximations. We neglect nuclear cor- 
rections to free structure functions in deuterium [§]. We neglect heavy-quark (charm and 
heavier) contributions. We assume qsea{x) = (jseaix). Finally, while displaying various equa- 
tions we ignore sea-sea contributions; these will be included in actual calculations, but are 
neglected here for simplicity. 

A. Drell-Yan Asymmetry for Proton Beams 

For proton projectile we have 

i7,7(Xi,X2) = I f4K(xi)«?(x2) +W?(X1)<(X2)) +<(X1)<(X2) + J?(X1)<(X2)1 , (2.4) 
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and a similar expression for ifi;^ (xi,X2). The asymmetry is 



A 



ttDY _ ttDY 

DY pp pn 



P RDY _|_ ffDY 

pp ' pn 



_ - (4<(xi) - rfg(xi)) (J^(X2) - u^,{x,)) + {Au^^jx,) - Jg(Xi)) «(X2) - dPjx^)) 
(4mS(xi) + rfS(xi)) [d^six2) + u^six2)) + (4uf (xi) + rf§(a;i)) (nS(x2) + rfS(x2)) 

In Fig. |1], we show the asymmetry yl^^(xi = X2) vs. xi. Recall that xi = Xheam and X2 = 
Xtarget- The curves are obtained with the phenomenological nucleon parton distributions 
of Gliick, Reya and Vogt [Q with Q^ = 20 GeV^. (Calculations were also performed with 
CTEQ-3M |l^ and MRS(G) |[Tl| parton distributions; results were similar and are not 



shown.) The lower curve is obtained with the full parton distributions, while the upper one 
results from fixing uP{x) = d^{x). For this kinematics, the two asymmetries have opposite 
signs: If we allow flavor asymmetry, the predicted DY asymmetry is negative, while it is 
positive if we require SU(2) flavor symmetry. In some kinematic regions the asymmetry Af*^ 
is large, but the DY cross sections may be so small that it cannot be accurately measured. 
Figure |l| requires that we measure the DY asymmetry along the line xi = X2- We can 
examine whether this is the most sensitive region to test flavor asymmetry in the nucleon 
sea. In Fig. ^ we show how the asymmetry A^^ varies with x\ and x^- In Fig. ^ja, we plot 
equal- Af'^ contours. Here and in other figures, the two extreme contours are labelled; the 
intermediate contours are understood to have labels differing in steps of 0.1. In Fig. [pD, we 
plot the equal- A^^ contours obtained by fixing viF{x) = dPi^x). In Fig. ||c we show contours 
of equal asymmetry difference Af*^ defined by 

Af^(a;i,X2) = <^|.=,-- <^|.^,- . (2.6) 

Along the line xi = X2 (the dashed curve in Fig. 2a), we see that the asymmetry A^"^ is 
negative and becomes more negative with increasing x. Secondly, if we set u^{x) = dP{x) 
and Xi = X2, the resulting asymmetry is positive and roughly constant at A^^ ^ 0.1 — 0.2; 
see Fig. 2b. Finally, the difference between the two asymmetries goes on increasing as 
Xi = X2 increases; see Fig. 2c. All these features are consistent with Fig. 1. Thus the 
equal-x kinematics is quite adequate if the purpose is to maximize the difference Af*^. 

However, if the purpose is to maximize the asymmetry A^^ , the equal-x kinematics does 
not seem to he the optimum choice. The asymmetry can he significantly altered hy deviating 
from the dotted line in Fig. 2a. If Xi = X2 = 0.2, the asymmetry is about —10% which is 
consistent with the experimental observation in Ref. 0]. However, if the measurements had 
been made at xi = 0.3 and X2 = 0.2, the asymmetry would have been ^ —20%. On the 
other hand, if the measurements had been made at Xi = 0.1 and X2 = 0.2, the asymmetry 
would have been nearly zero. Experimentally the ratio (o'bd/'^o'bp) is measured, and then 
the asymmetry A^^ is deduced. Contours for a fixed {(Jbd/'^o-bp) are identical to those for 
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a fixed Ap (Fig. 2a); only the contours labels are different. 

The E866 group at FNAL ||^ has measured DY processes for proton projectiles on protons 
and deuterons over a wide kinematic region. They focus especially on the region of large xi = 
Xbeam and small X2 = Xtarget- For Sufficiently large xi, the contribution from the projectile 
sea in Eq. ^^ is negligible relative to the valence contribution. Since d^{xi)/u^{xi) < 1, in 
this kinematic region the ratio of Drell-Yan cross sections is approximately given by 

as can be seen from Eq. ^^. The E866 group is thus able to measure the ratio d /vP over a 
wide range of x f^ . 

B. Drell-Yan Asymmetry for Kaon^ Beams 

For K~^ projectile we have 

H^Xp = ^ [4 {u^ ul + uT K) + dT dl + ^f " sl\ , (2.8) 

and a similar expression for H^Xn- Here we have suppressed for convenience the dependence 
of the parton distributions on Xi or X2, since there is no ambiguity in this case. The 
asymmetry is 

ttDY _ tjDY 
.DY _ ^K+p ^K+n 

""" ~ H^X. + H^X. 



K+p "I" ^-"K+n 



(2.9) 



4nf + (Jf + M?) + (4Mf + + d^^) {uZ + dl) + 2sf + si ' 

In Fig. ^ we show the asymmetry A^{xi = X2) vs. xi. The curves use the GRV 
parton distributions for the nucleons 0. Since little experimental information is available 
on kaon structure functions, we use in their place the GRV distributions for the pion 0. 
The lower curve uses the full parton distributions, while the upper curve is obtained if we 
fix u^{x) = d^{x). Note that for x between 0.1 and 0.4, the difference between the two 



asymmetries is quite large. If one had an intense well separated K^ beam, it would, in 
principle, be an excellent probe of flavor asymmetry in the proton sea. 

As is clear from Eq. |2]^, in order to calculate A^ we need four different parton distri- 
bution functions for K^ , while the only experimental information we have is on the ratio 
u^ /u^ W^. This ratio was found to be consistent with unity for xi < 0.7, and was less 



than unity only for larger values of xi. It was determined by using an unseparated negative 
beam containing pions, kaons and antiprotons, and by making the following assumptions: 
the meson sea was neglected, among the valence quarks in the mesons only the up flavor was 
retained, and finally the Drell-Yan scale factors for K~ and 7r~ were assumed to be iden- 
tical and constant over the kinematic range explored. This points out the urgent need for 
more experimental data on kaon structure functions, and justifies our use of pion structure 
functions in place of kaon structure functions in the calculation of the Drell-Yan asymmetry. 

We again examine whether the region of equal-x kinematics is the most sensitive region 
to probe the flavor asymmetry in the proton sea. Figures 4 a-c correspond to Figs. 2 a-c, 
respectively, with the projectile proton replaced by a K^ projectile. In the case of the proton 
projectile, we saw that the equal-x kinematics was adequate for the purpose of maximizing 
the difference Af*^ (Fig- 2c), but was not optimum for the purpose of maximizing the 
asymmetry A^"^ (Fig- 2a). In the case of the K^ projectile, we find that the equal-x 
kinematics is not good for either purpose (Figs. 4a and 4c). It is clear from Fig. 4a that if 
Xi = X2 < 0.3 the asymmetry is nearly zero, but can be significantly enhanced by deviating 
from the Xi = X2 line. 

The asymmetry for K~ mesons incident on protons and neutrons (in the deuterium) can 
be calculated similarly. Using the same approximations as before, we obtain 



ttDY _ 1 r 
^K-p - g 



4 (uf " < + ur < + ur <) + dr < + sr -si\ , (2.10) 



and a similar expression for Hj^-^. The asymmetry is 

ttDY _ ttDY 
aDY ^ "-K-p ^K-n 
K- ttDY 1 ztDY 



-4m 



K- 



dl-vFA + Uu'} +4Mf -d 



-,K- 



jK- 



« - dl) 



f2.111 



4u^- (£ + ul) + (4Mf - + 4Mf- + Jf -) (mS + rfS) + 2sf-s?' 
Numerical calculations show that this asymmetry is insensitive to the up/down difference 
in the proton sea-quarlc distributions (see section IV). 



C. Drell-Yan Asymmetry for Pion Beams 



For 7r+ projectile we have 



H. 



DY 

n+p 



1 r 
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4 {uf ul + «f ul) + cif dl + Jf rf^ + df dl 



and a similar expression for H^/^. The asymmetry is 





ttDY 

aDY -"tt+p 


ttDY 


dl) 




■' V+ rrDy _|_ ttDY 

^TT+p + -"TT + n 

_ - {Auf - Jf ) (rf> - u^) + (4tzf - df - df) {ul - 


IT 


{Auf 
beams, one obtc 

ttDY 1 


+ df) [£ + uf) + [Auf + df + df) {ul + 
dns 

4 («r w^ + wf ul + uf ul) + df dl + df dl 


il) 

1 



and a similar expression for Hj^f. The asymmetry is 



DY 



H. 



DY 

■K~p 



H. 



DY 



ffDY + jjDY 

■K p ' TT n 

Aul- - dt) (dl - uf + (Au:- + Auf - dl 



\u'. 



dl) 



(2.12) 



(2.13) 



(2.14) 



A 



{Aul + dl ) (dl + ul) + [Aul + 4mJ + rfj ) {ul + dl) 

Numerical calculations show that both A^J and A^Y are insensitive to the up/down 
ence in the proton sea-quark distributions (see section IV). 



(2.15) 
differ- 



Ill. DRELL-YAN SUM RULES 



If one assumes (a) charge conjugation symmetry so that e.g., q^ = ^ , (b) isospin 
symmetry so that e.g., u"" = d'^ , (c) qsea{x) = <isea{x)-, and (d) that charm and heavier 



flavors are negligible, then it is straightforward to show that tt"*" has only three charged- 
parton distribution functions (PDF), namely 

Uval{x) = d^al{x) = V (x) , 

"•seayXj UggnyX) "seayXj "seayXj = ^UghtyXj, 

Ssea{x) = Ssea{x) = ^ heavy \^)-i (3-1) 

and those for tt^ can be identified with one of these three. Thus vr^ and tt^ together have 
only three charged-parton distribution functions. With the same set of assumptions (a)- 
(d), it is easy to show that K^ ^ K^ , K^ and K^ together have only five charged-parton 
distribution functions which we denote by 

^light\^)i ^heavy\^)i ^light,l\^)i ^light,2\^)i ^heavy\^)i (3-2) 

where the notation is as in Eq. (3.1). 

To experimentally determine the three PDF's for pions, Eq. (3.1), one needs data on 
any three of the following four quantities: 

ttDY ttDY ttDY ttDY /o o\ 

^n-Dy -"yr-p! -"vr+D? ^n+p- l-J'-jJ 

Out of these, the first three may be preferable due to higher event rates. In fact, if the 



sea-sea contribution to H^"^ is ignored then Hj^"^ would contain only V^i^x) and S^^, 



light 



X 



and data on only two of the four quantities in (3.3) would be sufficient. 



Similarly, if the sea-sea contribution to H^^ is ignored, then Sf^^^^y{x) does not occur in 
H^^ . To determine the remaining four PDF's for K, one needs data on 

ttDY ttDY ttDY ttDY (oa\ 

^K-D^ ^K-pi ^K+Di ^K+p- l-J-^J 

If further the difference between Slfgf^^^ and 5*^^^^ 2 is ignored, then the data on any three 
quantities in (3.4) would be sufficient. 

As an alternative, some of these PDF's can be determined with the help of DY sum 
rules which we now present. For a beam particle B incident on the deuterium (D) target, 
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we assume H^}^ = H^^ + H^^ . Recall also that xi = Xbeam and X2 = x tar get- It is 
straightforward to show that 



iH!^~l - 


ttDY\ 


dX2 


= 7/9.y-(xi), 


(i^.% - 


ttDY \ 


dX2 


= V-{x,), 


(H^'^P - 


- ^ir,) 


dX2 


= 8/9 . \4j,,(xi) 



J {H^'^j.-H^Xn) rfx2=4/3.\4j,,(xi), (3.5) 

where the limits of integration are and 1. Note the interesting fact that these rules do 
not require any knowledge of the target particle structure functions. If the integration was 
done over xi instead of X2, we would have got a different set of sum rules with the well 
known nucleon structure functions, instead of the poorly known meson structure functions, 
occurring on the right-hand sides of Eqs. (p.5|) . 

There are additional sum rules where the integration is done over both x\ and X2'- 

J J {Hg^~H^^) dxi dx2 = l7/9, 

{H^^~H^^) dxi dx2 = 3, 

{H^'^j, - H^Xp) dxi dx2 = 8/9, 

(H^'^o-HJ^Xd) dxi dx2=A/3, 

{H^y^-H%) dx, dx2 = 7/9, 

{H^^^-H^^^) rfxi dx2 = l. (3.6) 

If QCD corrections are included, all the right-hand sides in the above sum rules will get 
multiplied by the well known DY K factor whose value is about 2. Since the nucleon PDF's 
are rather well known, one may also look at the first of Eqs. ( p.6| ) as a way to determine 
the average DY K factor. It may not be easy to verify these sum rules experimentally, 
because each rule involves two experiments (e.g., with vr" and vr"'' beams). Moreover, the 
integration is over the entire kinematic range which would entail some extrapolation of the 
data. However, if the beams do not have great intensity, the integration may help in getting 

10 



enough counts to make a measurement. Thus, if suitable beams are available, experimental 
verification may not be impossible. Sum rules similar to those in Eq. (|3.6|) were given by 



Hwa |T3|. The fact that the valence contributions to the DY process can be isolated by 
taking the difference between particle beam and antiparticle beam on a given target was 
mentioned by Sarma WM, but no sum rule was given. 



IV. CONCLUSIONS 

Why are the proton and K^ projectiles more sensitive than the K~ , vr"*" and 7t~ projec- 
tiles, to the flavor asymmetry in the proton sea? Inspection of the expressions for A^^ for 
various projectiles shows that they all contain a contribution proportional to {d^{x2) —u^{x2)) 
and a contribution proportional to {u^{x2) — d^{x2))- The former constitutes the "signal" 
while the latter is "noise". For all five projectiles, the signal is getting multiplied by valence 
parton distribution functions in the projectile. But these projectiles differ from each other 
in the following respect. In the case of p and K~^ projectiles, the noise is getting multiplied 
by the sea parton distribution functions from the projectile, while for the K~ , Tr"*" and vr" 
projectiles, the noise is getting multiplied by valence parton distribution functions from the 
projectile. Naturally in the latter case, the signal is expected to be relatively weak. This is 
borne out by the numerical calculations presented here. The difference between the numer- 
ical results for p and K^ projectiles on one hand and the K^ , 7r+ and tt^ projectiles on the 
other, can be traced to the fact that p and K~^ contain neither u nor d as valence partons, 
while K~ , 7r+ and tt^ do. 

Results presented in Figs. 2 and 4 bring into focus the kinematic regions which can 
be optimally explored in future experiments designed to probe the flavor asymmetry in the 
proton sea. Finally we have presented sets of Drell-Yan sum rules; those in Eq. (3.5), for 
example, could throw some light on the poorly known meson structure functions. The first 
of Eqs. (3.6) could be used to deduce the DY K factor averaged over the full kinematic 
region. 
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FIGURES 
FIG. 1. Theoretical Drell-Yan asymmetry Ap {xi = X2) vs. xi for proton beam; see Eq. (2.5). 
Here xi = Xbeam and X2 = Xtarget- The lower curve is calculated with full parton distributions, the 
upper one fixing vP{x) = dP{x). 

FIG. 2. Contour plots of Drell-Yan asymmetries for the proton beam, (a) Contours of equal 
asymmetry Ap , using full GRV parton distributions of Ref. [9]. (b) Contours of equal asym- 
metry using GRV parton distributions but fixing vP{x) = dP{x). (c) Contours of equal difference 
A^^, defined in Eq. (2.6). The two extreme contours are labelled; the intermediate contours are 
understood to have labels differing in steps of 0.1. Dashed curve: Xbeam = xtarget- 

FIG. 3. Same as Fig. 1 except that the proton beam is replaced by the K^ beam. 

FIG. 4. Same as Fig. 2 except that the proton beam is replaced by the K^ beam. 
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